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Abstract 



The reformulation of nonequilibirum thermodynamics, to include the treatment of thermody- 
namic fluctuations, is applied to the hydrodynamic fluctuations of a simple fluid. It is shown 
that the nonequilibrium thermodynamic scheme leads to the explicit form of the Fokker-Planck 
equation which describes the time behaviour of the probability distribution function of these 
hydrodynamic fluctuations as well as the irreversible processes which are connected with this 
behaviour. 



1 Introduction 



In a recent paper [1] one of the authors (P.M.) showed that the scheme and formahsm of nonequi- 
hbirium thermodynamics can be reformulated in such a way that the equations describing the 
dynamics of thermodynamic fluctuations are obtained in a simple manner. In particular one is 
thus led directly to a theory of fluctuations as Gaussian Markov processes obeying a multivariate 
Fokker-Planck equation [2,3]. In this reformulation of thermodynamics of irreversible processes 
fluctuations were considered as internal degrees of freedom of a thermodynamic system. More- 
over Gibbs entropy postulate was used to define the system's entropy as a functional of the 
probability density in internal coordinate space. 

In this paper we apply the formalism developed in [1] to the case of hydrodynamic fluctua- 
tions in a simple fluid. It is our aim to show that nonequilibrium thermodynamics yields indeed 
the multivariate Fokker-Planck equation describing the behaviour in time of hydrodynamic fluc- 
tuations and thus implicitely also the Langevin-like equations of fluctuating hydrodynamics [4]. 

In Section 2 we introduce the hydrodynamic fluctuations as fluctuation of the Fourier trans- 
formed (in space) flelds of mass, density, temperature and the three components of fluid velocity. 
We also introduce in this section the probabihty density function of these fluctuations as well as 
the system's entropy as a functional of this distribution function. This functional for the entropy 
which is an expression for Gibb's entropy postulate contains also the equilibrium distribution 
function. We give an explicite calculation in Gaussian approximation of the latter in section 3 
(and the appendix). In Section 4 we calculate the entropy production when the system is not in 
equilibrium. We obtain for this quantity an expression which is, in analogy to its form in macro- 
scopic nonequilibrium thermodynamic [5], a sum of products of dissipative fluxes, which are 
components of the probability flux, and conjugate thermodynamic forces, which are derivatives 
of probability densities. As usual in thermodynamics of irreversible processes linear relations 
containing phenomenological coefficients are established between fluxes and forces. When the- 
ses are introduced into the continuity equation for the time dependent probability density the 
differential equation, the Fokker- Planck equation, governing the dynamics of the hydrodynamic 
fluctuations as Markov processes, is obtained. We then determine in section 5 the phenomeno- 
logical, or Onsager coefficients, introduced formally in section 4 by deriving mean regression laws 
from the Fokker-Planck equation which was found, and comparing these with the macroscopic 
equations for the corresponding irreversible processes. 

We conclude with a few remarks, made also previously regarding the applicability of the 
formalism used. 

2 Hydrodynamic fluctuations as internal degrees of free- 
dom of a thermodynamic system 

We consider a simple fluid enclosed in a cubic container of volume V, and having constant 
total mass M and total energy E. The macroscopic thermodynamic internal state at time t 
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and position r, of the system may be characterized by 5 fields Ai{f,t) corresponding to the 5 
conserved quantities p(f, t), e(r, t), pv{r, t) where p is the mass density, e the energy per unit of 
volume and v the local velocity. 

For our purpose it is more convenient to use an alternative and equivalent set of fields which 
we choose as follows 

A,if,t) = pif,t) 
A2{f,t) = T{f,t) 

Ai{f,t) = v„(f,i), (i = 3,4,5; a^x,y,z), (2.1) 

where T(r, t) is the temperature field and where the index a denotes the 3 cartesian components 
of V. In the equilibrium state the fields A^ have uniform constant values A^ 

Al^Po = y; Al^To; A = 0, i = 3,4,5 (2.2) 

It is around these values that the fields Ai{f,t) will fiuctuate due to thermal agitation. To 
deal with these fluctuations a statistical description is needed, a statistical description which can 
best be effectuated in term of the Fourier components of the fields (2.1). To introduce these we 
choose for simplicity's sake periodic boundary conditions and write 

A{f,t)^V-'^e''^-%{k,t) (2.3) 
it 

with 

A{k,t)= / dre*A(r,t) (2.4) 

The wave vector k in eq. (2.3) and (2.4) will assume a denumerable set of values. 

We shall now introduce the probability density P{{Ai{k)},t) that the system is at time t 

— * 

in a state defined by the set of values {Ai(k)} = A. To normalize this probability density we 
introduce a cut-off wave number fcmax such that no state occurs for which \k\ > k^ax- The 
wave number k,nax in fact defines a length scale which is much larger than the molecular scale 
but small enough for a hydrodynamic and thermodynamic description to hold. With the above 
considerations in mind we now normalize the probability density as follows 

P{A,t)dA^l (2.5) 

I k I ^kmax 

In the limit as t — > oo the distribution function P{A,t) becomes the equilibrium distribution 
function Pq{A) for which the state variables have mean values {Ai{k)) 

(Aik)) = J A{k)PoiA)dA (2.6) 
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Substituting eq. (2.4) into eq. (2.6) we have 

{A{k)) = y(A(r1)%,o (2-7) 
where we have used the fact that {Ai{r)) is independent of r and that 

h,o-y'"jdr^'' (2-8) 

Since the equiUbrium averages of Aj(r) should be equal to the values of the equilibrium state 
given in eq. (2.2) it follows that 

(Ai(^))=M%o; {A,{k))^V%5^,,- (A,(^)) = 0, i = 3,4,5. (2.9) 

We now define fiuctuations 

ai{k) = Ai{k) - {Ai{k)) (2.10) 
These obey the distribution functions P{a,t) and Po{a,t) introduced above. We also note at 

— * — * 

this point that the fiuctuations ai{k) vanish for A; = 0: while this is true rigorously for ai as 
a consequence of mass conservation it is also true to linear order in fiuctuations (that is to the 
order relevant for our developments) as a consequence of energy and momentum conservation. 

We may now consider the fluctuations ai{k) as internal degrees of freedom of the fluid in the 
sense of nonequilibrium thermodynamics [6] as was done in the more formal preceding paper 
[1]. In this perspective in which the fluid is considered as a many component system with mass 
densities MP{a,t), its total entropy S obeys a Gibb's equation 

5S^-M [ da^5P + ^SE - ^5M (2.11) 
J To Tq Tq 

where /x(a) is the chemical potential of the component with mass density MP{a,t) and /^o the 
chemical potential per unit of mass of the fluid system at equilibrium. 

The entropy S{t) can moreover be expressed as a functional of P{a,t) according to Gibb's 
entropy postulate 

r P(a t) 

S^-kB P{a, t) In „V ^ da + So{E, M) (2.12) 

J Po[a) 

where 5*0 is the entropy at equilibrium and ks Boltzmann's constant. This expression yields for 
the differential of S 

SS = -kB I cia In ^^^5P + ^5E - ^SM (2.13) 

Comparing eq. (2.11) and eq. (2.13) and using eq. (2.5) we flnd for /x(a) 

/ V kBTn , P(a, t) , ,s 
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The previous formulae in particular cqs. (2.11) and (2.13) on which our further analysis will be 
based indicate the essential role played by the equilibrium distribution function Po{a). It is for 
this reason that we shall now first establish the explicit form of this function for the variables 

— * 

ai{k) introduced above. 



3 The equilibrium distribution function of fluctuations 

To discuss the equilibrium distribution function Po{a) wc use the connexion between the prob- 
ability of a given state a and the value of its entropy S{a) 

Po{a) = const x e^^^^/^^ (3.1) 
Po{a) = Po(0)e^^(^)/*^« (3.2) 



or 



with 

AS{a) = S{a) - S{0) (3.3) 

where S{0) refers to the entropy of the equilibrium state. 

We will therefore now calculate AS{a) to second order in the fluctuating variables which we 
have introduced in section 2. We first observe that the entropy of our fluid may be written as 
the volume integral of an entropy density per unit of volume s(r) 



AS = f As{r)df (3.4) 



At each point r the entropy density s{r) and its fluctuation A will be a function of the 
internal energy density u{r) and the mass density p{r). 

If we then calculate A^* to second order in Am and Ap which is sufficient within the context 
of a linear fluctuation theory we have 

AuAp+lf^(ApA (3.5) 



dudp 2 dp^ J 

In the appendix, we show that this expression can be reduced to a quadratic form in the 
fluctuations Ap(r), AT{r) and v{f) or alternatively Ap{k), AT{k) and v{k) 

AS = -4/(if{^(|)To(Ap)2 + ^{ATf + pov^} 
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^AT(fc) AT(-fc) + pov{k) ■ v{-k)} (3.6) 
where p is the pressure and c„ = pQ^du/dT the specific heat per unit of mass; P — {dP/dp)T 

— * 

denotes the isothermal compressibility. The fluctuating variables in the sum over k correspond 
to the five variables ai{k), i = 1,2, ■■■ ,5 introduced in section 2. 

We shall redefine in a minor way three of these variables, namely those for i = 3,4,5 and 
rewrite the scalar product v{k) ■ v{—k) as follows 

v{k) ■ v{-k) = v{k)v{-k) + iF{k) ■ if{-k) (3.7) 

where the scalar function v(k) and the transverse vector 'iF(k) are given by 

v{k) = k ■ v{k) (3.8) 

{/^(fc) = (1 - kk) ■ v{k) (3.9) 
with k the unit vector in the direction of k, k = k/k. Defining now the modified variables 

a^{k)^v{k) (3.10) 



a,(A;) =<(A;), (z = 4, 5; 7 = 1, 2), (3.11) 
with 7 denoting the transverse orthogonal components of v^{k), can be written as 

= -\j:g,,{k)a,{k)a,{-k) (3.12) 

k,ij 

The real matrix gij{k) in this quadratic form is diagonal, gij — guSij and has elements 

9u = {VTopo)-'l3; ^22 = {VT^y'poC, 
gn^{VTo)-^Po; i = 3,4,5 



(3.13) 



Formulae (3.2), (3.12) and (3.13) completely determine the equilibrium distribution function 
Po{a) and yield for variables Xi{k), conjugate to a.t{k) according to 

X,{k) ^ -kB^^^ = j:9iMk), (3.14) 
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t liG rclcXt iQxis 

X,{k) = {ToVpo)-'PAp{k), = {VTi)-^poc,AT{k) ] 

(3.15) 

Xi{k)^{ToV)-^Pov{k), i = 3,4,5 J 

We shall now evaluate in the next section the entropy produced in the system when the 
hydrodynamic fluctuations do not have their equilibrium distribution and derive, following the 
rules of non-equilibrium thermodynamics, the Fokker-Planck equation describing the dynamics 
of these fluctuations. 



4 Entropy production and the Fokker-Planck equation 
for hydrodynamic fluctuations 

Since the system considered is a closed one with constant total mass and energy, its rate of change 
of entropy with respect to time is according to equations (2.11) and (2.13), or equivalently (2.12) 

— = -kB daln \ / — 4.1 
at J Po{q.) ot 

We shall now make use of the continuity equation obeyed by the probability density P(a, t) 

^-^--H^Ju^PM (4.2) 

— * 

where J^^ can be interpreted as a current per unit of phase space density along the k,i axis 
of the phase space spanned by the complex state variables aj:-. Equation (4.2) represents as it 
were, the local form in a-space of probability conservation. Notice that, since Ai{r) is a real 
quantity, the state variables ai{k) obey the relation a*{k) = ai{—k). Consequently, volume 
elements in the phase-space will be real and positive and since probabilities are real and positive 
by deflnition so will be the probability density P{a) as it should. Equation (4.1) then implies 
that the currents Jk,i have the property Ji. = 

Wc now substitute eq. (4.2) into (4.1) and obtain, integrating by parts and using the fact 
that P{a,t) tends to zero sufficiently rapidly for ai{k) tending to infinity, for the irreversible 
rate of change of entropy 

dS r 

— = P{a, t)a{a, t)da > (4.3) 

with 
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the local entropy production in a-space. As usual in non-equilibrium thermodynamics the en- 
tropy production is a sum of products of currents and thermodynamic forces. In the present 
case the latter are the gradients in a-space of the chemical potential (2.14). 

The total entropy production for the system is according to eq. (4.3) an average over possible 
fluctuating states of a (a), the entropy production accompanying the evolution of a particular 
state. Since moreover the inequality (4.3) is valid for an arbitrary choice of the probability 
density P{a), it follows that the second law of thermodynamics holds not only for the average 
of a (a) but also for a (a) itself 

(7(a) > (4.5) 

In equilibrium when a{a) vanishes, both the currents and thermodynamics forces occurring in 
expression (4.4) vanish. In non-equilibrium and within a sufficiently wide range of physical 
conditions we may then suppose that these quantities are connected by linear relations. In 
a strictly linear theory of hydrodynamic fluctuations to which we will restrict ourselves, the 
coefficients in these laws, linear in the thermodynamic forces, must be constant with respect to 
the variables ai{k), while the various modes do not couple. The linear laws will then be of the 
form 

— * — * 

If we furthermore introduce the variables Xj{k) conjugate to aj{k) according to eq. (3.14) 
relations (4.6) can alternatively be written as 

Jk,i = - E Lij (Xj + kB^^^ In P(a)) (4.7) 

Substituting equation (4.7) into the conservation law (4.2) the following differential equation 
is obtained 

which describes the time evolution of the probability density P{a, t). 

With the linear relation between Xi{k) and aj{k) (c.f. eq.(3.14)), the differential equation 
(4.8) can also be written in the form 

where 

= J2 LiiQij, = Yl ^iiafj' (4-10) 
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and 

Lg = ^(L,, + L,,) (4.11) 

This equation has the standard form of a hnear Fokker-Planck equation and describes there- 
fore the dynamics of the hydrodynamic fluctuations as Gaussian- Markov processes. As stated 
previously this result has been obtained here within the framework of non-equilibrium thermo- 
dynamics and its extension to system with internal degrees of freedom. 

5 Identification of the Onsager coefficients in the Fokker- 
Planck equation 

Note that the Fokker-Planck equation (4.9) contains the scheme of coefficients Mjj or equiva- 
lently Lij (the matrix gij connecting these has already been specificied above) which must still 
be determined to make the description of the dynamics of the fluctuations complete. The de- 
termination can be achieved by observing that fluctuations regress in the mean according to the 
macroscopic laws for irreversible processes. This experimentally verified bahaviour is known as 
Onsager 's regression hypothesis [7] and enables one to identify, in the case of linear laws the 
coefficients Lij in terms of macroscopic quantities. 

We establish the mean regression laws from the Fokker-Planck equation (4.9) by multiplying 
both of its members by ai{k) and integrating over all variables a. We then find for the first 
moments of the variables ctj after partial integration the ordinary differential equations 



^^^^ = -Y.L^Ak)xr\Kt) 

= -Y.Mu{k)a^\k,t), (5.1) 
I 

where the conditional mean is given by 

a^'' ^ J ai{k,t)P{a,t)da (5.2) 

with 

P{a, 0) = n <^ - «io(^)) (5.3) 

k,i 

The mean regression equations constitute one of the fundaments of Onsager's derivation [7] 
of reciprocal relations for the coefficients L^j, which follow by applying microscopic reversibility. 
For the equations (5.1) containing complex variables ai{k), defined through equations (2.1) these 
reciprocal relation are [8] for complex Onsager coefficients L^j, 

Lij = L^iSiSj (5.4) 
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with 

\ -1 for i = 3,4,5 ^^-^^ 

The coefficients satisfy, in addition to the reciprocal relations, certain conditions imposed 
by the spatial symmetry of the system. For the isotropic fluid considered this implies that scalar 
and vectorial phenomena do not couple so that 

L,, = L,, = for{^.^ ^1^,3 ^^ g^ 

and also that the matrix L^j, i, j = 4, 5 is diagonal 

Lij^LSij, ij^4:,5 (5.7) 
To identify finally the coefficients we hst below the famihar [9] fully linear macroscopic 

— * — * 

hydrodynamic equations for the density, temperature and velocity fields, p{k,t), T{k,t) and 

— * 

v{k,t). These equations are 

dp{k,t) 



dt 



— —ipokk ■ v{k,t) (5-8) 



mM = .^iML,a . ,Ck, t) - ^r(S, t) (5.9) 

dt poCy poCy 

dv{k,t) ik(3 - ik ( dp\ rj^if 

r.. = - P{k,t) [^\ T{k,t) 

dt po po \dTJp^ 

- y^^k^v+[lv + Vv^kk-]v{k,t) (5.10) 

In equation (5.9) A is the heat conductivity; r) and r)v in equation (5.10) are the viscosity and 
volume viscosity respectively. 

By multiplying moreover eq. (5.10) succesively with the unit vector k, and the projection 
operator 1 — kk one obtains the macroscopic equations for the variables v{k) and v^{k) introduced 
in section 3, cf. eqs. (3.8) and (3.9). 

Comparing then the mean regression equation (5.1) with the macroscopic hydrodynamic 
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equations one finds for tlie matrix of phenomenological coefficients M, 



^3 



M 









ikpo 





















POCv 


pocv yar) 


ik§_ 
po 


ik ( dP\ 
PO \dT)po 



















k^ri 
po 

















k^rj 



(5.11) 



PO 



and for tlie matrix of Onsager coefficients L^, using also eqs. (3.13) and (4.10), 
/ ikVTo 



VT,'k'X/plcl ik{%)pVTi/plc, 
tkVTo ik{%)pVTi/plc, VTok\\r] + r],)lpl 












V 



















VTo^rj/pl 













(5.12) 



VTok'rj/pl J 



Tfiis tfien concludes the specification within the context of nonequilibrium thermodynamics 
of the Onsager coefficients Lij for the linear fluctuating fluid considered and completely deter- 
mines the Gaussian solution of the linear Fokkcr-Planck equation. Observe that the reciprocal 
relations (5.4) as well as the spatial symmetry requirements (Curie's principle) (5.6) and (5.7) 
are automatically satisfied by the coefficients found above. 

The above application to hydrodynamic fluctuations of the recently developed extension 
to fluctuation phenomena [1] of the scheme of macroscopic nonequilibrium thermodynamics, 
illustrates once more previously stated characteristics of this generalized formalism which goes 
beyond the traditional macroscopic domain. It introduces indeed, while keeping the framework 
and rules of the macroscopic theory, statistical notions such as: 

1. The distribution function of hydrodynamic fluctuations. 

2. The definition of entropy as a functional of this distribution function. 

3. The calculation of the probability distribution function in equilibrium. 

We also stressed before that the extension of the thermodynamic formalism is legitimate 
because it is limited to events taking place on a time scale which is slow compared to the time 
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scale of microscopic phenomena. For the case at hand the phenomena dealt with occur at 
characteristic times which are essentially hydrodynamic relaxation times. 



Appendix 

In this appendix we show for completeness sake that expression (3.5) for AS, the variation 
of entropy with respect to its equilibrium value up to second order, reduces to the diagonal 
quadratic form in AT(r) and Ap(r) given in formula (3.6). 

We first observe that in the integral (3.5) the linear terms give rise, due to mass and energy 
conservation to the integral of the quadratic term ^pov'^{r): 

since all thermodynamic derivatives in (A.l) are uniform equilibrium quantities and since ds/du — 
I/Tq. Next we observe that the second order terms in the integral (3.5) may be written in the 
following form 

as 
and 

T dp dpdu dp^ 



-A^ = = + ^ Ap. {AA) 



Expand now in expression (A. 2) Am and Ap/T to linear order in AT and Ap 

du , ^ du ^ , . . 

^^-^^^^0-p^P (^-^^ 

A^ = ^Ar + l|^Ap (A6) 

T dT Tdp ^ ^ ' 

Here again the thermodynamic derivatives are uniform equilibrium quantities. 

Finally observe that the total differential of the thermodynamic potential //T, with / = 

u — Ts the free energy density, 

d^^ = -ZidT+^dp (A7) 



leads to the Maxwell relation 

d/i/T 1 du 



{A. 



dT T2 dp 

Substituting then (A.l) and (A. 2) with (A. 5) and (A. 6) into expression (3.5) and taking into 
account the Maxwell relation (A. 8) one obtains the desired result (3.6). 
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